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Abstract 



We propose and analyze an improved small-x equation which incorporates exact 
leading and next-to-leading BFKL kernels on one hand and renormalization group 
constraints in the relevant collinear limits on the other. We work out in detail the 
recently proposed w-expansion of the solution, derive the Green's function factor- 
ization properties and discuss both the gluon anomalous dimension and the hard 
pomeron. The resummed results are stable, nearly renormalization-scheme inde- 
pendent, and join smoothly with the fixed order perturbative regime. Two critical 
hard pomeron exponents uJciQ"^) and uJsiQ'^) are provided, which — for reasonable 
strong-coupling extrapolations — are argued to provide bounds on the pomeron 
intercept wp. 
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1 Introduction 



Recent results on the next-to-leading log s corrections |T|,|^ to the BFKL equation and 
to the hard pomeron show subleading effects which are so large as to question the very 
meaning of the high-energy expansion and thus raise the compelling question of how to 
improve it. 

Two facts suggest that an essential ingredient of any improvement of the BFKL ap- 
proach should be the correct treatment of the collinear behavior, as predicted by the 
renormalization group (RG): on the one hand the success of normal QCD evolution in 
explaining the dependence of the small-x behavior of structure functions at HERA, 
and on the other hand the observation that the large next-to-leading corrections to the 
BFKL equation come mostly from coUinearly-enhanced physical contributions. 

A first attempt at introducing collinear improvements was performed long ago, by 
the treatment of coherence effects 0] in the collinear region. This leads to the CCFM 
equation |^,^], which differs from the BFKL equation by subleading effects to all orders, 
even if a full inclusion of the DGLAP splitting functions in a consistent CCFM framework 
has not yet been achieved. Other modified BFKL approaches incorporating some DGLAP 
evolution are being tried too 

Very recently, it has been suggested |^ that such all-order collinear effects can be 
incorporated as subleading kernels of a generalized equation, whose solution can be found 
by the method of the so called cj-expansion, allowing in particular a resummation of the 
energy-scale-dependent terms of the kernel . 

The purpose of the present paper is to insert such suggestions in a general scheme, 
which leads to the renormalization group improved small-x equation, and to study its 
solutions. One of the outcomes will be to stabilize, in a nearly scheme-independent way, 
the estimate of the anomalous dimensions and of the Q^-dependent hard pomeron. 

A first point to clear up is which pomeron we are going to estimate. Previous work on 
RG factorization |10|,|11|] in the BFKL equation with running coupling |r^-[T9| has shown 



that the pomeron c<jp defined as the Q^-independent leading singularity in the u;-plane, is 
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beyond the reach of a strictly perturbative approach. On the other hand, there appears 
to be a boundary of vahdity of the RG, the Q^-dependent hard pomeron uj^lQ"^), which 
is argued to be independent of the small-fc^ strong couphng region and is thus hopefully 
calculable in perturbative QCD. 

Since ujp{Q'^) signals a change of asymptotic regime, it is associated with an u- 
singularity of the anomalous dimensions, not necessarily of the full gluon distribution. 
Thus LJp{Q'^) may be related to a power-like behavior in an intermediate small-x, moder- 
ate region, and not to the very small-x asymptotic behavior of the structure functions. 
It also follows that uj]f>{Q'^) is a rather difficult quantity to determine, because it is related 
to the position of an cj-singularity, and is thus dependent on the full anomalous dimen- 
sion perturbative series. Possible definitions, leading to a precise estimate, are thoroughly 
discussed in Sec. |^ and in Sec. H, where our results are provided. 

A second point to realize is that, in order to incorporate the collinear behavior cor- 
rectly, a whole string of subleading kernels, represented by a series in the running cou- 
pling as(Q^), is to be taken into account. In fact, the leading logs calculations count 
one high-energy gluon exchange per power of ag, with any transverse momentum ratios. 
In the collinear limit, provided by the strong ordering in the transverse momenta, only 
the singular part ~ 1/z of the DGLAP splitting function is obtained. The remaining 
part contributes to higher and higher order subleading kernels which carry fewer powers 
of logs, but leading powers of logQ^. From a quantitative point of view, such collinear 
contributions are very important, and in fact account for most of the exact NL kernel 
(Sec. 3). 

Let us stress that we are not aiming here at a full control of subleading log s contri- 
butions, but only of those that carry a leading collinear contribution from NNL level on. 
Therefore, we remain in a context in which only t-channel iteration is important, without 
mixing with the s-channel one (see, e.g., [0,^). 

In this framework, we can define the fc-dependent gluon distribution by the NL k- 
factorization formula introduced by one of us in large k dijet production in parton- 
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parton scattering. In a general hard process involving probes A and B we can write 
daAB [ du f s 



dfc dfc, 







/ ^ (^) ^^(^)(^l(l + "si^) 6^.. {l + a,H^)\ko)hB{ko), (1.1) 



= Ifcl , /Cn = Ifcfl 



where the impact factors h^, hs may carry additional dependence on the hard scales of 
the probes and the gluon Green's function is provided by 

g^{k,k^) = {k\[uj - Kj\-^\k,) , (1.2) 

apart from the multiplicative kernels H, W which may be needed at subleading level O, 
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We notice immediately that the scale of the energy s in Eq. ( |1 . 1| ) has been chosen to 
be kkf), i.e., factorized and symmetrical in the "upper" (u) scale k and "lower" (/) scale 
kQ. This means that Q^u and the kernel /C^; in Eq- (|1-2| ) are both symmetrical operators. 
On the other hand, when k ^ ko {ko ^ k), the variable kko/s is not the correct scaling 
variable, which is rather k'^/s {kQ/s) — i.e. the usual Bjorken variable. 

In order to switch to, say, the upper energy-scale fc^, it is apparent from Eq. ( |1 . 1| ) that 
one has to perform a similarity transformation {k/k^j^Q^, which in turn implies 

the relationship 

lC^:\k,k') (^^J = ICUk,k')=IC^\k,k') (^^y (1.3) 

between the symmetrical kernel /Co; and the kernel /C^"-* (/C^^-*). Although technical, this 
remark is important in order to classify the collinear logarithms, because if a wrong 
energy-scale is chosen, single logs (of k/k') may turn into double logs (cf. Sec. 3). 

The main purpose of our study is to construct the RG improved kernel /C^, and to 

provide the solution for Q^^ in Eq. ( |1.2| ) in the RG regime k^ ^ k"^ ^ A^. The starting 

point is the observation p| that the kernel }C^{k, k'; fi'^; as(/i^)), for non-vanishing values 

of k, k', is RG invariant, and can thus be expanded as a power series in as(fc^) with scale 

invariant coefficients 

oo AT ^ 

/c.(fc, k') = j:Me)r' K:ik, fc') , = ^ = ^1^^^^ • (1-4) 
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Since we want to take into account the leading collinear singularities to all-orders, the 
series ( p. .41) is necessarily infinite, as noticed before. 

Solving for the Green's function with the general kernel ( |1.4|) is a novel problem in 
the BFKL approach, which is addressed and solved in Sec. 0. There we derive the main 
properties of the solutions, namely (i) the factorization property of Q^j in the RG regime, 
(ii) the u- expansion of the relevant eigenfunctions and {iii) the definitions of the pomeron 
singularity tup and of the hard pomeron singularity ti;p(as(^^))- In a first reading of this 
rather mathematical section one could perhaps retain the basic results, and come back to 
their derivations after Sees. 3-5. 

In Sec. 3 we explicitly construct the improved kernel /C^^ with the requirements of [i) 
reducing to the exact L-I-NL terms in the relevant limit and {ii) reproducing the known 
collinear singularities at higher orders. 

The corresponding solution for in the RG regime and the explicit form of the solu- 
tion JF^ of the homogeneous equation are studied in Sec. 4. The main result is that the NL 
truncation of the improved tu-expansion takes into account correctly all collinear singu- 
larities, at least for the purely gluonic case. The inclusion of the (small) qq contributions 
is discussed also. 

Finally, in Sec. |^ we present our results for the resummed anomalous dimensions 
and for the hard pomeron, and exhibit their stability under scheme change and NNL 
corrections. 

In the last Sec. ^ we discuss the present situation and future prospects, which in- 
clude a solvable model |^ , based on the collinear analysis of the present paper. A few 



mathematical details are covered in the Appendix. 



2 Small-x equation for a general kernel 

We consider here a general form of the small-x equation, whose as(t)-dependence is sup- 
posed to be consistent with leading-twist anomalous dimensions and must contain, there- 
fore, an infinite series of subleading terms (cf. Introduction). Our final goal is to investigate 
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the solution for the gluon Green's function (|1.2| ), i.e., the resolvent of the improved kernel 
LuG^ik, ko) = 6\k -ko)+ /C^(fc, k')g^ik', ko) , (2.1) 

in order to derive its large-t behavior in the RG regime. 

2.1 Form of the kernel 

The improved kernel /C(^(fc, k') occurring in Eq. is assumed to have the asymptotic 
as (t) -expansion 

}CUk,k') = J2[asit)r^' Kik,k') , t = log^ (2.2) 

n=0 

where the coefficient kernels are scale- invariant and may be cj-dependent. They are 
partly known in closed form from leading and next-to-leading [|I],0 calculations, and 
have known ^ collinear properties to all orders. 

The leading coefficient kernel must reduce, for — 0, to the historical one 
having eigenvalue function 

Xo(7) = 2^(1) - ^ii) - V^(l - 7) , ^ = Y^ (2.3) 

on test functions {k'^y^^. The NL coefficient kernel K'^ is related also, for u ^ 0, to the 
one recently found on the basis of NL QCD vertices, except for the subtraction of a 
term already included in the ^-dependence of Kq (cf. Sec. 3). 

In general, the expansion ( p.2|) was justified in Ref. |^ as follows. First IC^{k,k'), 
at energy-scale sq = kko (Eq. ( |1 . 1| ) ) and non vanishing virtualities, is a collinear finite 
distribution, symmetrical in its arguments. By RG equations, for k and k' much larger 
than A, JC^^ must have the form 

/C,(fe,fc';/i2;a3(/i')) = ^/C^(«s(t); t, = ^/C^(«s(t'); i', i) (2-4) 
which, by expanding in as(t), yields Eq. (|2.2| ). 
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In the limit of vanishing virtuaUties (A; ^ or A;' ^ 0) IC^) acquires colhnear singu- 
larities, which are dictated by the nonsingular part of the gluon anomalous dimension in 
the Qo-scheme which, by neglecting the qq part, is 

7(dj) = -igg{uj) - — = a,Ai{uo) + cilA2{uj) + ■ ■ • , (2.5) 
Ai{uj) = + 0{u) , A2{u) = Q + 0{uj) . 

As a consequence, the eigenvalue functions Xni.l) acquire the 7-singularities 
X-(7)^ ^-^-^^-t^^,-V^f.Vr^^-^^^ (7«1) 



(7 + 

l-{Ai-h){Ai-2b)---{Ai-nh) 



1 - 7 « 1) , (2.6) 



(1 - 7 + i^)"+' 
where h is the one-loop beta function coefficient (cf. Sec. 3). 

The 7, uj dependences are tied up together in Eq. (p.6|) because of the similarity 
relations ( |1.3| ), which define the kernels /C''"-* (/C*^'^) at energy-scale k"^ (/cq) having simple 
collinear behavior for k ^ k' {k' k). As a consequence, the 7-singularities occur at 
shifted values of 7 (by ±co'/2) and the symmetry of /C^^ imphes, by Eq. ( p. 21) , a slightly 
asymmetrical 6-dependence in Eq. ( p.6|) . 



2.2 Factorization of non-perturbative effects 



In order to actually solve Eq. (2J.) for ^^(fc,fco), one should extend the representation 
( ]2.2| ) in the region around the Landau pole A; ~ A (t = 0), where it becomes unreliable. 
Whether such an extension can be somehow hinted at on perturbative grounds — as 



in the time-like evolution case p4| — is an open problem that we do not address here. 



However, for perturbation theory to be applicable, the non-perturbative effects of such 
region should be factorized out, as is predicted by the RG, and has been argued for at 
leading and first subleading level [pl],pT 



In the following, we consider the dependence of Q^^ on various kinds of regularization 
of K,uj in Eq- (P-2|) around the Landau pole, and we argue that indeed the RG factorization 
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property holds for sufficiently large t, in the form 



Q^{k, ko) = J^^{k)J^^{ko) + higher twist terms , (t — to > 1) • (2.7) 

Here JF^ (JF^) is the solution of the homogeneous small-x equation 

/C^^c. = ujj^u. , (2.8) 

which is "regular" for t +oo [t —>■ — oo) in the sense that it is asymptotically in the 
corresponding region (see Sec. ^]3|for a more precise discussion). 



Let us ffist try to understand how Eq. (|2.7] ) can possibly work. By inserting it in the 
defining equations 

uj9u.it, to) - j dt' xut, t') g^(t', to) = S{t - to) , 

X^{t,t') = kk'IC^{k,k'), SUt,to) = kkog^{k,k'), (2.9a) 

and by using the symmetry of Sw, we obtain, for t — to ^ 1, 



uS'^it) - I dt' XUt,t') 3^^{t') ~ (2.9b) 

A^(t',to) 



OO 

to 



dt'X^{t,t') 



oo 



3^^{t') 5^ - ^^(f) 



-Fujito) 

?^(t) = fc^^(fc), 3^^{t) = k:F^{k 



dt' X^{t,t' 



3^Lu{to) 



where A^(t, to) denotes the higher twist part of in Eq. (|2.7| ). Now, let us go to the 
large-t limit: the l.h.s. is the homogeneous small-x equation for 3^^^, and the r.h.s. will be 
negligible, i.e., higher twist, by the following mechanism. First, note that t' < O (to) in 
the r.h.s., because A^, by definition, decreases rapidly for |t' — to| ^ 1. Furthermore, for 
|t — t'l ^ 1, JCuj satisfies the coUinear factorization of Sec. 3 (with higher twist corrections), 
so that the to-dependence in the r.h.s. is factored out and can be made to vanish by a 
proper choice of A^. 

We thus conclude that, provided the regularization of the running coupling allows 



such properties of the kernel, the factorization in Eq. ( p.7| ) of the large-t dependence 



actually holds. The decomposition of the kernel in a factorizable and in a local part is 



certainly satisfied in the case of models leading to differential equations (cf. Ref. JTTl and 



the collinear model of Ref. as soluble examples), but is presumably satisfied also in 
the case of kernels in an C? space having reasonable spectral properties, as we shall argue 
next. 



2.3 Form of the solution 

We thus assume that, by a suitable regularization of aj^^ around the Landau pole, K.^ 
can be defined as a hermitian operator bounded from above in an C? Hilbert space, with 
a continuum (or possibly discrete) spectrum — oo < /i < ii^ioS). Typical regularizations 
of this kind may (a) cut-off a^if) below some value t = t > 0, or (b) freeze it in the form 
aJ(C) = (6t)^^0(t — t) + (6f)~^0(t — t), possibly with some smoothing out around the cusp. 
The spectrum of /C^^ is expected to be discrete in case (a) and continuum in case (b) |Tl|] . 
In the latter case, the expansion in Eq. ( |2.2| ), extended to the region t < i, defines a scale- 
invariant kernel with frozen coupling, where however the coefficient kernels K!^ should 
be evaluated, for consistency, in the 6 = limit. This limit introduces some ambiguity 
in the definition of below i, which in our point of view is part of the regularization 
procedure. 

In such a framework, a formal solution for the Green's function is given by the 
spectral representation 

in terms of the eigenfunctions 

/C^J-^ = /i^^, (2.11) 
which satisfy an orthonormality condition 

(-^.^-^.^') = / ^ JT(fc)^.^'(fc) = Jdt J^^*m^:{t) = 6{f, - /i') (2.12) 
and can be chosen to be real, because so is /Co; (A;, k'). 
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We shall normally consider the situation for which Re(a;) > fiw{uj), so that uj is not a 
point of the spectrum (|2.11|) , and JF^ (JF^) in Eqs. ( |2.7|) and ( |2.8| ) are not eigenf unctions, 
being well behaved for t ^ +00 (t — — 00) only. 

We shall also refer, in most of this section, to the example of the frozen-Os regular- 
ization, which allows a simple classification of the eigenfunctions J^l^{k) of Eq. 
according to their behavior for t ^ —00. In fact, since the test functions 

= (^2^-igi.(M)i ^ (1 _ ^ = 1/2 + iz/) , (2.13) 



are reproduced for large negative t by the kernel ( |2.2| ) with eigenvalues 



/i = X][as(t)]"+^Xn(l/2 + H^). b = 0) (2.14) 



n=0 



the eigenfunctions Tl^{k) must have the behavior 

jr''(fc) = 1 (F'^(^)(fc) -F-'^('^)(fc)) (2.15) 
~ 1- (r(z/)e''^(^)* - r*(z/)e-'"(^)*) 

for suitable functions F^{k) having a plane- wave asymptotic behavior for large and nega- 
tive t (the uo index has been dropped). The two "frequencies" z/(/i) and —t^ifi) correspond 
to the two solutions of Eq. ( |2.14| ) for real fi, which are real also, because of the 7 1 — 7 
symmetry of Xnil) ^^e 6 = limit, as better seen from their explicit form, similar to 
the basic one in Eq. ( |2.3|) (Sec. 3). Note also that the spectrum endpoint is provided by 
the maximum of the (real) expression (|2.14D when u varies. 



The precise superposition of left- and right-moving waves occurring in Eq. ( |2.15| ) is 
determined by the condition that jF^(fc) be regular for t —>■ +00, i.e., be vanishing at least 
as rapidly as 1/k, so as to allow an (continuum) normalization. While the negative-t 
behavior ( |2.15| ) is oscillating for /i in the spectrum ( p.l4| ), it becomes a superposition of 
decreasing and increasing exponentials when /i is continued off the real axis with Re(iz/) > 



0. This structure, similar to that of potential scattering [11|, suggests that the Green's 



function can be asymptotically evaluated by the "on shell" expression 

GUk, fco) ^ -^"(fc) F^^^H^o) , t - to » 1 , (2.16) 
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thus identifying J^uj{ko) = F^^'^\ko) in Eq. ( |2.7] ) as the solution of the homogeneous BFKL 
equation which is regular for to ^ — oo. 

The argument goes as follows. By using Eq. (|2.15|) we rewrite the spectral represen- 
tation (|2.10| ) as a contour integral 

where F'^^'^'* and F"'^^^^ are assumed to be boundary values of an imaginary analytic 
function of fi, whose branch cut lies along the spectrum and is encircled by the contour 
C{uj). We then evaluate the behavior of ( p.l7[ ) for to — oo, by distorting the /x-contour 
(because F'^^^^ is well behaved, for Re(iz/) > 0) and by picking up the residue at the ^ = uj 
pole, i.e., the r.h.s. of Eq. ( |2.16|) . This procedure can be carried through for to > also, 
where F'^^^'^ becomes irregular, provided t — to is large enough for the decrease of to 
compensate the increase of F^'^^^ 

The plausibility argument above is further supported by the explicit model of Ref 



for arbitrary values of t and to, and hints at the general validity of Eq. ( |2.16|) . Therefore, 
for k ^ /cq, the Green's function is asymptotically proportional to the "on shell" regular 
solution of the homogeneous BFKL equation J-'^{k), which becomes the basic quantity 
to be found. 

Furthermore, the above procedure allows us to define also the pomeron singularity 
u = Up. In fact, the integral representation ( ^.17| ) is singular when its contour is pinched 
between the branch-point /i = /ip(u;) and the pole fJ^ = w, i.e., for i^(/ip = a;p) = 0, or 

oo 

up = /ip(^p) = ^[a,(tl]"+^xr(l/2, b = 0) (2.18) 

n=0 

which is an implicit equation for up in the present regularization procedure of ag-freezing 
at small k. For a general regularization, the definition up = fip{ujp) is still valid, but the 
explicit expression (|]T|) is not. 

It follows that ujp is a singularity of the right-moving wave F^^^^ rather than the regular 
solution, and that it affects the asymptotic behavior (|2.16|) in the to-dependent coefficient 
only. Therefore, the regularization dependence of up and of the spectrum is factorized 
away asymptotically. This picture is confirmed by the explicit examples of Refs. [0,PB 
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2.4 Small-cj expansion 



We follow the philosophy of Ref . , according to which <^ 1 is the relevant expansion 
parameter of the solution, rather than as{t). Furthermore, we first consider the "off-shell" 
case fi u, or more precisely fi < fip{uj) < Re(ci;) <^ 1, and we take the generalized ansatz 

(2.19) 

where X^(7,/i) is to be found by solving Eq. (p.ll|) . 



Once again, we are interested in the RG regime 6t > ^ 1, in which the regular 
solution in Eq. (|2.19| ) turns out to be dominated by the stable saddle point 7 = 7^(/i, t) 
defined by 



{7^-^^0.(7,/")} =0 



7=7 

-xL(7,/^)>0, (#)'^|-(#). (2.20) 

It has already been shown ^ that, around the saddle point, the effective eigenvalue 
function Xuj{l,fJ') is independent of the regularization procedure and its /x-expansion has 
been evaluated for /i = tu, by a treatment of the saddle point fluctuations (App. [A. 11 ). 



Here we prefer to find the /i-expansion, in the same regime, by using the replacement 
t ^ d.y in 7-space [p!0|p!^ , in order to give an all-order evaluation. We thus write Eq. ( |2.11| ) 



for t > t in the form 

1 .-i. 1 



f!^ = i^lCo + —fifC, + ^^IC, + ...jfli (2.21) 

and by repeated partial integrations we prove the 7-space identity 

&M(^?(7)/.^(7)) = [(X.(7,/") -M7)^?(7)]/.^(7). (2.22) 



Strictly speaking, the validity of Eq. (|2.22|) is limited by the fact that the bfit operator 



has to be regularized around t = (e.g., by freezing it for t < t). However, the large-t 
behavior of Eq. ( |2.19D can be safely evaluated by (|2.22|) provided 



fe/it ~Xc.(7,/i) (2.23) 
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by the saddle point condition ( p.20| ). 

By replacing ( |2.22| ) into ( p.21|) we obtain the equation 



X^(7,/^) = Xo(7) + (X" - b^^^,r ' ^^Xlil) + (x" - ^/^^t)" V'x2 (7) + 



(2.24) 



which, at a given subleading order in Xn provides a nonlinear differential equation for 
Xuj{l,f^), and thus a formal solution of Eq. ( ^^.211 ). 

However, since we are looking at the large-t and small-/i limits, we prefer to expand 
Eq. ( p. 24 ) in the denominators as well, thus obtaining the following asymptotic expansion: 



X.(7,/i) = Xo (7) + /^<(7) + /i'%"(7) + ■ 
where (App. [70| ) 



(2.25) 



V3 = —; 







xr 




1 


VXo 


x^ 


1 


' X^3 b 



xr 

x^ 



(2.26) 



X2 



xTxl 
(x^)^ 



I 7, o;/ n uj LO 



and so on. This expansion is supposed to yield safely the large-t behavior of Eq. ( p.l9| ), 
whenever Eq. ( p.23|) is satisfied. The cumbersome saddle point fluctuation method of 
App. |A.l checks with the result in Eq. ( |2.26|) . 



The /x-expansion of the regular solution in Eqs. ( p.25| , 2.26 ) is the basic result of this 
section and will be applied in the following ones to actual NL calculations. 



2.5 Anomalous dimension and hard pomeron 

Due to the validity of RG factorization in the large-t limit of Eq. ( |2.7| ), we can state that 
the gluon density g^{t) in the probe A has a universal t-dependence 

9^it) = ^M^) - k'g^ik, fco) ~ k'T.ik) , (2.27) 



13 



where the A-dependent coefficient is in general non-perturbative. By Eq. (p.l9|) , this 
yields the proportionality relation 



where we have specified the function at the "upper" energy-scale Sq = A;^ (7 ^ 
7 — ^uj), which is relevant in the large /c^ limit. 

The asymptotic behavior of Eq. (|2.28| ) in the RG regime can be found from the saddle 
point ( p.20| ), which yields the result 

^ ( / .A r-TT7T^ + ■ ■ ■ ) exp / 7(0;, r) dr , (2.29) 
where 7(0;, t) = %{uj,t) satisfies the identity 

7t - X^") (7, 00) = [ 7(cu, r) dr + const (2.30) 



and the coefficient in front, coming from the saddle point fluctuations, has been evaluated 
at NL level only. 



If we work at NL level, the saddle point approximation ( ^.291 ) is enough, and provides 



the effective anomalous dimension [25 



■rM^.t) = *) - ( i + 1^^^^ 1 + . . . (2.31) 

XW'(7,u;) \^7 2 x^"^'(7,^^) 



whose subleading expansion has however a wildly oscillating behavior The hard 

pomeron singularity comes in this case from the failure of the saddle point expansion at 
the point uj = ujs{t)), such that 

X {'y{uJs,t),LJs) = (saddle point estimate) , (2.32) 

thus implying infinite fiuctuations in Eq. (|2.31|) . 



On the other hand, in our RG improved approach, we do not rely on a subleading 
hierarchy. Therefore, the estimate ( p.32| ) may be not realistic. For instance, it has been 
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suggested that Eq. (|2.31|) yields higher-order singularities of oscillating type which 
may perhaps resum to a scale change. Here we just notice that Eq. ( |2.28|) , with the solution 



( |2.25| ) can simply be evaluated beyond the saddle point approximation for oj < Ug, and 
yields a generalized definition of the effective anomalous dimension 

7efr(^, t) = -^7- = . (2.33) 



The analysis of the ^-singularities of Eq. (|2.33|) has been done in some toy models 



in which the BFKL equation reduces to a differential one |TT|, ^]. In such cases the 
singularity comes just from a zero of the denominator at = uJc{t), so that 

9u)c.{i)if) = (gluon-density-zero estimate) . (2.34) 

The two definitions ( [^.32| ) and ( |2.34| ) yield in a sense two extreme estimates of the 
hard pomeron singularity {udt) < ujp(t) < uJs(t)) of the full anomalous dimension series. 
Both will be discussed here on the basis of our improved BFKL kernel. 



3 Improved subleading kernels 

The general form (|2.2| ) of the kernel of the RG improved small-x equation is strongly 
constrained by (i) the exact leading and next-to-leading logs calculations and (ii) 
the collinear singularity structure of Eq. ( p.6|) . This leads to a natural identification of the 
coefficient kernels and — up to some NNL ambiguity — following the procedure 
of Refs. which is described in detail here. 

3.1 Form of the collinear singularities 

Let us first recall the argument leading to Eq. ( |2.(j| ). The RG invariant kernel in Eq. ( |2.4D 
acquires collinear singularities in the limit k'/k — > {k/k' —>■ 0), which corresponds 
to strong ordering of the transverse momenta in the direction of the "upper" scale fc^ 
("lower" scale k^). Therefore, such singularities are easily expressed for the kernel /C*-"-* 
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(/C^'^) corresponding to energy-scale k"^ (fcg) in the NL fc-factorization formula (|1 . 1|) . For 
k ^ k', /C*^"-* acquires the form 



lC^:\k, k') ^ ^ exp 7(^, «s(r)) dr (t » t') 



p ^l-&«s(t)log^J ^ , (3.1) 

where 7 is the non-singular part of the gluon anomalous dimension of Eq. ( |2.5| ), the 
singular one being accounted for by the BFKL iteration itself. 

Expanding Eq. ( p.l|) in as{t) and comparing with the general definition ( p.2|) , leads to 
the identification of the kernels K^"*^ in the collinear limit, whose eigenvalue functions 
turn out to have the singularities 

+ + (3.2) 

which correspond to single logarithmic scaling violations for k ^ k^. A similar reasoning 
yields the collinear behavior of /C'^'-''^ in the opposite strong ordering region k' k 



K^\k,k') c^^expj^ 7(^,«s(r)) dr {f » t) 



"^^^^ ^-6a.(t)log^l ' (3.3) 



k'^ \ ' ^ k^ 

and to the singularities 

However, the similarity relation ( |1.3| ) connects the kernels /C*^"-* and /C^^-*. Therefore 
/C(") has the singularities (|3.4| ) shifted at 7 = 1 + to also, and similarly K,^^^ has the 
singularities (|3.2|) shifted at 7 = — cu. As a consequence, the symmetrical kernel K.^ — for 
the energy-scale Sq = kk^ — has both kinds of singularities shifted by ±co'/2, as anticipated 
in Eq. (|2.6| ). In particular the leading and NL coefficient kernels have singularities 



X^(7) - + z \^ (3.5) 

(7+2^) (1-7+2^) 
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Note the 6-dependent asymmetry of the singularities in Eq. ( |2.6|) under the 7^1 — 7 
transformation. It is due to the fact that the expansion (|2.2|) involves as{t) (and not 
as(t')). Of course, the kernel itself must be symmetrical under t ^ t' exchange, so 
that expressing dsit') in terms of asit) 

= "^^'^ (3.7) 

l + bas{t) log-p 

leads to the symmetry constraints 

^"(^) = E (^) i-bd.r-'^X^n.a - 7) , (3.8) 
where denotes the 7-derivative. It is straightforward to check by the binomial identity 

that the symmetry constraints ( p.8|) are indeed satisfied by Eq. ( p.6| ). In particular we 
must have 

Xj^(l - 7) = Xo (7) , xr(l - 7) = xr(7) + bxo'h) , (3.10) 
showing that the antisymmetric part of Xiil) is ~|xo'(7)- 

3.2 Form of the leading coefficient kernel 

Given the fact that the w-dependence is tied up with the 7-dependence in the singularities 



( |2.6| ), it follows that the leading logs hierarchy, corresponding to a pure tu-expansion at 
fixed 7, is poorly convergent close to 7 = and 7 = 1. This observation follows from the 
trivial expansion 

' ifl-^+(^V+--| (3.11) 



7 + |cj 7 \ 27 V27 / 

and was used in Ref. to suggest a resummed form of the leading kernel eigenvalue 
function 

X^(7) = [^(1) - ^(7 + M + [^(1) - ^(1 - 7 + Ico)] (3.12) 

=xo{v - 7;^— — + ••• • 

2 sm 7r7 
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The kernel Kq , corresponding to Eq. ( p.l2D is that occurring in the Lund model and 
is given by 



(3.13) 



where fc> = Max {k, k') and k^ = Min {k, k'). It is thus related to the customary leading 
kernel Kq by the "threshold factor" (k^/ky)'^. This means that the s-dependence provided 
by its inverse Mellin transform is 



Kois; k, k')^ I (^)" K-ik, k') = Koik, k') e(s - k 



(3.14) 



Can one justify the form of the kernel (|3.13| ) "a priori"? From the point of view of the 
RG improved equation, any kernel which (i) reduces to Kq in the — > limit and (ii) 
has the leading simple poles of Eq. (p.6|) for n = 0, is an acceptable starting point. An 



alternative choice of this kind will differ from Kq by a NL kernel without 7 = or 7 = 1 
singularities. The resulting ambiguity can thus be reabsorbed by a proper subtraction in 
the NL coefficient kernel. 



Nevertheless, the threshold interpretation of Eqs. ( p.l3| ) and ( p.l4| ) is appealing. For 
instance, the first iteration of such a kernel provides the expression 

/ , , N /" du; / s f 1 ^ ^ 
KQ{s;ki,k2) 



27ri \k1k2 
d^k 



UJ 



(3.15) 



Ko{k,,k){\og 



kik2 



r]{ki, k) - r]{k2, k) Ko{k, ^2 



where 



cosh?7(A;i, k) 



e + kf 



(3.16) 



The threshold condition implied by Eq. ( p.l5| ) 

g 

——— = cosh 77 > coshirjiki, k) + rjik, ^2)) 
2kik2 



(3.17) 



is reminiscent |30| of phase space in Toller variables |^ and may be regarded as an 
alternative way of stating coherence effects as implied in the original version of the 

Lund model itself . 
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Whether or not such hints will eventually provide a more direct justification of ^ 
the fact remains that Eq. ( |3.12|) resums the cu-dependence of the 7-singularities, and thus 
provides the correct singularities of the scale-dependent terms of the NL kernel. Therefore, 
it is a good starting point, yielding NL contributions which are smoother than those in 
the as(^)-expansion, as we now discuss. 



3.3 Form of the next-to-leading contribution 

The NL contribution K"^ is constructed by requiring that {€) the Green's function Q^^ 



reproduce the known NL calculations and {ii) the coUinear singularities be as in Eq. (2.6) 
with n = 1. 

In order to implement condition {i) we have first to relate the cu-dependent formulation 



of Qi^ in Eq.(L2) to the customary expression of the BFKL kernel at NL level 



-lC = -ia,Ko + alKi + ---) . (3.18) 



The cu-dependent formulation of Eq. (|2.2| ) yields instead the NL expansion 

i/C- =- (a,K^^^ + a.uK'i^ + alxf^ + ■■■). (3-19) 



Kf = Kf^ + uKI^^ + ■ 



which is actually more general than Eq. ( |3.18| ) because the a^oJ term, coming from the 
cu-expansion of Kq, is a possible NL contribution too. 



Now it turns out that, at NL level, the formulation ( p.l9|) reduces to the one in ( p.l8 



provided the impact-factor kernels if, H'^ of Eq. ( p. . 1| ) are taken into account. In fact, by 
using the expansion ( p.l9| ) and simple operator identities, we can write 

1 \ — i / 1 \ — i 



UJ J V / \ a; 



(3.20) 



provided we set 

= , ^1 = + I {<'Ko + Koir«) . (3.21) 
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Eqs. ( p. 201) and ( |1 . 1|) show that the two formulations above differ by just a redefinition 
of the impact-factor kernels, while Eq. ( |3.21|) means that Kf^ is given by i^i, after 
subtraction of the term already accounted for in the tu-dependence of Kq. Using Eq. (|3.12|) 
this yields the a; = limit of the eigenvalue function 



xr°(7) =Xi(7) + ^Xo(7)-4^ 
2 sm 7r7 



(3.22) 



The subtraction term so obtained is important because it has cubic poles at 7 = 0, 1 



which cancel the corresponding ones occurring in the energy-scale dependent part 



XoJ 



of Xi{l) found by Camici and one of us 0, as noticed by another one of us and 
seen explicitly in Eq. (|3.23|) . Furthermore, the impact-factor kernels of Eq. (|3.2CI| ) have 
quadratic poles which similarly account for the ones occurring in H and H'^ 0j0- This 
means that the remaining contributions are, in both cases, much smoother in the ui- 
dependent formulation. 



In order to implement condition (m) on Xii we note that the uj = Q limit ( |3.22D still 
contains double and single poles at 7 = 0, 1, which should be shifted according to Eq. ( p.6|) . 
By neglecting the (small) qq contributions, the explicit form of Eq. (|3.22| ), following from 
Refs for the energy-scale Sq = kkQ, is 



I (§(Xo(7) + Xo(7)) ) + 



-7X0(7) + 7^X0(7)—^ 

4 2 sm TT'j 



+ 



TT 



+ 



cos'7r7 
4 |Vsin-7r7y 3(1 — 27) 

67 7r2\ _ 3 



11 + 



7(1-7) 



l + 27)(3-27), 
3 



+ 



TT 



36 



1X0(7) + :^C(3) + -^ 

12/ 2 4sm7r7 



$(7) 



(3.23) 



$(7)=E( 



n=0 



V'(n + l + 7) - V(l) , ipin + 2 - -f) - 



{n + 7)2 



+ 



{n + 1- -i f 



Here we have singled out some singular terms which have a natural physical interpretation, 
namely the running coupling terms (in round brackets), the energy-scale-dependent terms 
(in square brackets) and the coUinear terms (in curly brackets). 

The running coupling terms have a double pole at 7 = 1 only, and account for the 
asymmetric part of xi (given by — |xo) which provides the 6-dependent double pole on Xi 
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in Eq. (|2.6| ). The coUinear terms have symmetrical double poles with residue Ai(iu = 0), 
in accordance with Eq. ( |2.6| ) also. Both types of singularities can be shifted by adding a 
NNL term, vanishing in the u = limit, which we take to be 

AH^'(7 + f ) - A(0)^'(7) + (A(^) - b)ij'{l - 7 + I) - (^(0) - b)ij'{l - 7) . 

(3.24) 

This term incorporates the u- dependence of the one-loop anomalous dimension ( p.5| ) too. 



The energy-scale-dependent term in square brackets contains the subtraction (|3.22|) 
and has, therefore, simple poles at 7 = 0, 1 only, which we can shift by adding the 
contribution 

y(Xo(7)-Xo(7)) . (3.25) 



By then collecting Eqs. ( |3.22| ), ( |3.24| ) and ( |3.25| ) we obtain the final eigenvalue function 

6 



xUl) = Xi(7) + A(^)^'(7 + i^) + (^i(^) - b)ij'il - 7 + 1^) + -Xo (7) , (3.26) 



XI (7) = XI (7) + 7;Xo{l)—T^ - -Xoil) - Ai(0)^'(7) - (^(0) - b)^\l - 7) 



where 

Ixoil)-^ ^ 

2 sm 7r7 o 

(3.27) 

is a symmetrical function without 7 = or 7 = 1 singularities at all. The expression 
( |3.26| ) satisfies in addition the symmetry constraints (|3.10|) , having antisymmetric part 



2X0 ■ 



Of course, there is some ambiguity involved in the choice of the subtraction terms 
(|3.24| , p.25|) , which boils down to the possibility of adding to ( p.27| ) a term, vanishing in 



the u = limit, and having only higher twist 7-singularities, around 7 = —1, —2, ■ ■ ■ and 
7 = 2, 3, ■ ■ ■ . This ambiguity leads to an error which is of the same order as that made 



in the NL truncation of the cu-expansion of the solution in Sec. |2.4| , as we shall see next 



3.4 Numerical importance of collinear effects at NLO 



Above we have given the general form for the collinear singularities of the kernel at all 
orders. It is of interest to consider at NLO just how much of the full corrections come from 
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Figure 1: A comparison of the coUinearly-enhanced (double and triple poles only) part of 
the NLO corrections with the full NLO corrections; nj = 0. 



these coUinearly enhanced terms. Accordingly we look at the part of the NLO corrections 
which contains just double and triple poles, xi,c- 

_ A,-b 1 

7^ + (1-7)2 273 2(1-7)3- ^ ' 

This is compared with the full Xi in figure where we have plotted their ratios to Xo- The 
remarkable observation is that over a range of 7, the coUinear approximation reproduces 
the true corrections to within 7%. It is obviously impossible to say whether this is true 
at higher orders as well. However the fact that the study of collinear terms has such 
predictive power at NLO is a non-trivial point in favour of our resummation approach. 



4 Improved next-to-leading solution 

Having constructed the coefficient kernels Kq and K'^ with consistent collinear behavior 
(cf. Eqs. ( ^.12| ) and ( pj.26D ) we would like to know the large-t behavior of the solutions of 
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the improved small-x equation, whose kernel (|2.2|) is truncated at NL leveL This problem 
has been solved in general in Sec. 2, and we describe here the NL features. 



4.1 cj- Expansion of the gluon distribution 

According to Sec. 2.4, the eigenf unctions J^^ik) (— oo < jj, < Hfiyj)) can be found in the 
small- /i, large-t regime 

fet > - > - » 1 (4.1) 



by the 7-representation (|2.19|) , i.e.. 



k'mk) = /'^'" :^exp Lt - ^XUl,A , (4.2a) 
ii-ioo 27ri bfi J 

where the exponent function is provided by the small-/i expansion 

d.XUl, /x) = xUl, /i) = Xo (7) + 4rT + ^2 (7) + /i' r/^(7) + • • • (4.2b) 

Xo wJ 

and , ''73 , ... are given in Eq. ( [^.2(j| ). 

Furthermore, by the factorization property ( |2.7] ), valid for |t — to| ^ 1? the gluon 
Green's function ( p.l| ) is itself proportional to Tuj{k) = T^^'^{k), which is obtained by 
setting II = uj m. Eqs. ( [4.2|) , i.e.. 



9,X^(7, u) ^ x(7, ^) = Xo (7) + ^ + NNL , (4.3) 

where we have now truncated the expansion to NL level. The ensuing error is argued to 
be small (Sec. [4.3|) . The RG regime holds if there is a stable saddle point 

but = x(7, ^ Xo (7) + ^ + NNL , x'il, u;) < , (4.4) 

Xo(7) 

which dominates the large-t behavior of Eq. ([4.2a|) , providing the anomalous dimension 
representation (|2.29 ). The effective anomalous dimension can be continued past its saddle 



point value by means of Eqs. ( |2.33| ) and ( p.28| ), which use the 7-representation ( |4.2a| ) for 
/i = a;. 
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4.2 Properties of the kernel and its solutions 



In this section we illustrate some of the features of the resummed kernel and of the regular 
solution as obtained with the 7-representation. 

It can be instructive to examine the resummed eigenvalue function xil^ ^) Eq. ([4.3|) 
in two different ways. Firstly as a function of 7 for various values of u, as shown in Fig. |^. 
In the symmetric energy-scale case, we can see how the original 7-poles at 7 = 0, 1 are 



so=kk' so = k^ 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 1.2 

r r 

Figure 2: ^(7, uj) as a function of 7 for various values of uj, for the symmetric energy-scale 
So = kk' on the left and for sq = k"^ on the right. Here n/ = 0. 



displaced in a symmetric way for 7^ 0. The slight asymmetry of ^(7, u;) is due to 
Eq. ( p.lOP , which in turn comes from the perturbative expansion in as{t), as discussed in 
Sec. Even for sizeable values of uj the eigenvalue function preserves its shape, with a 
stable minimum (x"(7m,C(j) > 0) around 7 = 1/2. This stability is a necessary condition 



to avoid the oscillating behavior noticed in Ref. 



With the "upper" scale choice sq = A;^, the pole at 7 = is tu-independent as it should, 
while the pole at 7 = 1 is shifted for 7^ 0. In this case we also have a stable minimum, 
in a slightly different position with respect to the previous case (7^ — > 7m + ■ 
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A second way of looking at the resummed kernel is to examine a quantity which we 
call xS(7,as), defined by 

xS(7,«s)=X^"H7,^ = «sxS). (4.5) 



This is closely related to the saddle-point approximation for evaluating the 7-representa- 
tion, ( ^.281 ), since the value of the saddle-point, 7 satisfies uj = asxig (7; ^^s); 7 such 
a case is itself closely related to the effective anomalous dimension, ( p.31|) . In Fig. ^ we 
show xl^il^Cis) for different values of as- The marked asymmetry is due to the energy- 
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Figure 3: xig (7!Q^s) as a function of 7 for various values of as, for energy-scale Sq = k'^ 
and Uf = 0. 



scale choice Sq = k"^- We note the rather different structure from the x(7,c<j) shown in 
Fig. ^ In particular there are no longer any divergences. That on the right is shifted by 
an amount uj: as a result rather than a pole one has Xeg ~ l/o^s, as discussed in 0. That 
on the left instead becomes Xcg ~ ase~"'^°'^ for negative 7 as a result of the inclusion of 
the dependence on the DGLAP splitting function (in particular the 1/(1 — z) part, which 
gives Ai{uj) ~ — logcu for uj +00). Another feature of xi!s "worth noting (though not 
immediately visible from Fig. ^ is that for n/ = 0, asxig (O^f^s) ~ independently of 
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as- This is so because close to 7 = 



(u) 
efT 



a. 



1 + ujAUuj) 



7 



(4.6) 



For 7 = 0, we have that 1 + ujAi{uj) = 0. Since v4i(l) = —1, = 1 at 7 = 0. This 
point is of significance insofar as it relates to the problem of ensuring conservation of the 
momentum sum-rule for the gluon distribution (its u = 1 moment). 

1 




Figure 4: gui{'t)e for u = 0.15, with energy-scale sq = fc^. The normalization is 
arbitrary. 



Next we examine the properties of the regular solution of the small-x equation for 



the integrated gluon distribution in Eq. (|2.28|) . In order to concentrate on the region 
most relevant for a consideration of the small-x properties of the anomalous dimension, 
in Fig. ^we actually show 5'(^(t)e~*/^, as a function of t = In A;^/A^. There are two critical 
points on the curve. Firstly the point labelled x'il) = 0, namely where the saddle-point 
solution but = x{l,^) sits at the minimum of x fo^' that u. Let us refer to that point 
as tg. The second point of interest, which we call tc, is the rightmost zero of g^j- This is 
the point where the effective anomalous dimension has a divergence. Since the solution 
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in this region has roughly the form [|13| 



2 



9^(t)-Ai\[-^] (bujt-x^m)]. (4.7) 
One can estimate the difference between tc and tg as being 

A.'n 



where — — —2.3381 is the position of the rightmost zero of the Airy function. For fixed 
as, this translates to a difference between the estimates for Ug and Uc which, for ^ 1, 
has the following form: 



(1,2 2 // \ J^/'J 
-^1^ j + O ^11.16 a^/^ . (4.9) 



Such af''^ contributions to oo have already been observed in other contexts where there is 
some form of cutoff on transverse momenta, such as a running coupling which is zero below 
a certain value of t, or non-forward elastic scattering (where the exchanged transverse 



momentum places an effective cutoff on transverse momenta in the evolution) ||l^|35 



Numerical estimates (based on the 7-representation) for the difference between Ug and Uc 
coincide with (f4.9|) but only for very small as- 

— —5/3 

For typical values of ag, we note that ujs ~ — 11-16q;s is of the same order as the 



NLL corrections. It too, as pointed out in |T^,|3^, is the first term of a poorly convergent 
series. The resummation procedure that we recommend (and adopt) is to define Uc not 
through the power series in as, but by looking for the rightmost zero of the regular 
solution. 



4.3 Estimate of error 

The question now arises: what is the error that we make in the NL truncation of the 
RG improved equation? Our claim is that, in the improved formulation, based on the 
cu-expansion ( [4. 3D , this error is smaller than in the formal NL expansion in as{t). Let us in 
fact estimate the remaining terms in Eq. ( [4.2b| ). According to Eq. (|2.6| ) further subleading 
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eigenvalue functions contain at least higher order collinear poles which contribute to 7/2 , 
?73 and so on. A first observation is that, even if has in + l)-th order poles, the 77^'s 
have at most simple poles, due to the powers of Xo in the denominator, roughly due to 
the replacement as{t) ~ ^/Xo- Therefore, their contribution cannot be too big, even for 
small values of 7 = O^u). 

Furthermore, one can check that, if qq contributions (Sec. 4.2) are neglected, the 



leading collinear poles actually cancel out in the expansions ( |2.26| ) of 773 , Vs^ • • • around 
both 7 = and 7 = 1. The mechanism of this cancellation can be cleared up as follows. 

From the mathematical point of view , it is possible to have the truncated NL solution 
to be an exact solution of Eq. ( p.8| ), provided the following recurrence relations hold 
(App. |A3) 

X2 f Xl Xl XS I Xl , rj \ X2 



^ - 69J ^ , ^ = l^-bd^]^, ■■■ . (4.10) 
Xo \Xo 7 X^ Xo \Xo V Xo 

It is now really simple to check that such relations build up the collinear singularities 
( |2.6| ), which therefore must cancel out in the subleading corrections r/g 5 ^73 5 • • • • The 
recurrence relations ( [4.10| ) can also be interpreted as DGLAP equations in 7-space, for 



the anomalous dimension 7 in Eq. (|2.5| ). 

From a more physical point of view, it is not possible for simple poles to survive in 7/2 , 
?73 , . . . because, when replaced in the saddle point condition (^.4|), they would provide 
cu^, uj^, . . . corrections to the one-loop anomalous dimensions which cannot possibly be 
there. In fact, the full anomalous dimension is accounted for by Eqs. ( [4.3| , [4.4|) as follows 

~ - + ^ = aj- + A.iuj)) , (4.11) 

7 7 / 

where we have taken the small-7 limit of the collinear safe eigenvalue function x^^\l,^)- 
We therefore conclude that, in the purely gluonic case, the NL cj-expansion (^.41) 



takes into account the collinear behavior to all-orders, and that no further resummation 
is needed. This point is perhaps more easily seen by replacing the NL truncation of 
Eq. ( [4.2b| ) in the saddle point condition ( p^.20[ ) to yield the equation 



b^it = x-o+^^^ /i = \. ■ (4.12) 
^0 1 - «s— 

x^ 
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It is apparent from the last version of Eq. ( [4. 121) that we are deahng with an effective 



eigenvalue function which resums the coUinear behavior as a geometric series. 



We are finally able to state that the error in the NL truncation ( 4.2b|) is uniformly 
0(0;^), the neglected coefficient having no 7 = nor 7 = 1 singularities at all. This error 
is therefore of the same size as the ambiguity in the definition of Xi f hat we have pointed 
out before. The corresponding error in the saddle point condition ( [4. 4] ) is a roughly 
7-independent change of scale A{bt) = 0{uj), or A{as) = 0{u)a^. 



4.4 Extension to qq contributions 

The coefficient kernels K!^ take up collinear singularities not only from the nonsingular 
part of the gluon anomalous dimension 7^^, but also from qq states which are coupled to 
it in the one-loop gluon/quark-sea anomalous dimension matrix 

Jabi^) = aAabiuj) = 7afe(t^) " ^ag'^^^ , (4.13) 

where a = {q,g) and Ca = {Cf,Ca) = Nc{r,l) denote the partonic channels and color 
charges. 

Although the numerical effect of quark-sea contributions to the gluon anomalous di- 
mensions is pretty small including the two-channel evolution ( [4. 13] ) changes the 



collinear problem conceptually. While the small-x equation stays of one-channel type, 
due to the high-energy gluon exchange, the two-channel collinear behavior yields two 
anomalous dimension eigenvalues 



7± = ± y (^^^V^ J +799799, (4.14) 

with the approximate NL expansions ([7sg/7gg] Leading 



7+ =^ 799 + ^"799 , 7- ^ 799 - ^^799 , ^ ^ = ) ' ^^'^^^ 



A 



Recovering in the BFKL framework the full collinear behavior ( ^.14| ) is not trivial, 
because 7_ ~ —r'jq^J = 0{as) starts at NL level and for 7 = 0{as) the leading logs 
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hierarchy breaks down in the ag-expansion What do things look hke in the uj- 



expansion? 

Note first that the derivation of the colhnear behavior of in Sec. 3 can be repeated, 
by replacing Ai{uj) with the matrix A(ct;) in Eq. ( ^4.13|) , and by projecting the final results 
onto the gluon channel, which corresponds to a bracket between initial state (^) and final 
state (r 1), because the quark couples to the high-energy gluon with relative strength 
r = Cf/Ca- Therefore, Eq. ( |2.6| ) should be replaced by 

Xnil) ^ ^l^y+A r l)|l.A.(A + 6).--(A + (n-l)&)|(^°^), 7«1, 

(4.16) 

^ n li .„^, ((rl)|l-(A-fe).--(A-n&)|f°')), 1-7«1. 

(1 - 7 + \ujY+' \IJ 



In particular 



(A) (A) - h 

(A(A + 6)) , ((A-6)(A-26)) 



Xl — , I 1 \2 ~^ /I I 1 \2 ' (4--'-''') 

(7 + ^^) (1-7 + ^w) 



where (A) = Agg + rAqg, (A^), . . . denote the brackets defined before in Eq. (|4.16| ). 



Secondly, the kernel ( |3.23D should be supplemented by the (gg) contribution , which 



completes the 6-factor in front of the running coupling terms and adds up a collinear 
contribution, as follows 



Xf(7)=-^(^(X^(7)+X[)(7))) + 



n 



f 



5 , , 3 TT^ cos7r7 l + |7(l-7) , 
Xo(7) + T7^ , ' L I ^ ■ (4.18) 



6Nc \3^^'^ sm'^TT-f 1 - 27 (l + 27)(3-27) 

Correspondingly, the subtraction term ( p.24| ) changes by the replacement^ 



Agg = ^ (A) = Agg + VA.g , (4.19) 

^Since there is a (small) two-loop anomalous dimension in the Qo-scheme, induced by qq contributions, 



one could envisage a shift of this simple pole in Eq. ( 4.18| ) also, by a further change of the NL subtraction 
term. 
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while Xo ^"^^ subtractions ( p.22D and (|3.25|) are left unchanged 



The main differences with the purely gluonic case come out in the cj-expansion of the 
solution, and specifically in the role of the higher-order terms. In fact, if we repeat the 
calculation ( [4.11|) with the new entries ( [4.17|) , we find 



tI'^ = + (A) j = «s + + rAggiuYj , (4.20) 

which is consistent with the NL expansion ( [4.15| ) for 7+, but is not the full one-loop 



anomalous dimension ( [4.14| ). 



Further terms in the tu-expansion must therefore contribute I/7 and 1/1 — 7 poles, 
and they indeed do. From Eqs. ( |2.31| ) and ( [4.17| ) we find 



A7« = ((A^) - (A)2) a^u + ((A - (A))^)^,.;^ + ■ ■ ■ (4.21) 



which checks with the explicit expansion of Eq. ( ^4.14|) up to the relevant order. The 



explicit matrix form of the corrections in Eq. ( [4. 21] ) makes it clear why the two-channel 



problem allows the survival of the simple 7-poles at higher subleading orders. 

Nevertheless, the small w-expansion remains smoother than the ag-expansion. In fact, 
the 0{uj^) NNL terms being neglected show simple poles only (around 7 = 0,1), the 
general trend remains the same as in Fig. |^, provided u is not too large. If uj increases, 7+ 
decreases, and at some critical value of uj, for which 7_|_ and 7_ become of the same order, 
the a;-expansion will break down, eventually. Whether or not the low-energy eigenvalue 
7_ can still be described by an all-order resummation in u remains an open question. 



5 Anomalous dimension and hard pomeron 

Here we present our main numerical results, for both the improved gluon anomalous 
dimension and the hard pomeron, and we show their stability. 
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Figure 5: The anomalous dimension in various approximations 



5.1 Results 

Fig. ^ shows the purely gluonic anomalous dimension as a function of u for cts = 0.2. 
The LL anomalous dimension is just •jll = Xo'^i^/^s) and has the familiar branch-cut at 
u = 4 In 2 The NLL anomalous dimension is taken as 

Xi{1ll) 



iNLL = 111 - ocs- 



(5.1) 



'Xo(7ll) 

and has the feature that it is always negative, with a divergent structure around the same 
point as 7^^. The resummed result, defined in Eqs. ( p^.28| ) and ( |2.33D , shows a divergence 



at a much lower uj, defined by uJcif) in Eq. ( 2.34 ). What is particularly remarkable is the 
similarity to the DGLAP result until very close to the divergence. The momentum sum 
rule is automatically conserved: for = 1 we have 7efT = (this is closely connected with 
the fact that cisxiff (0) — ^) — P^^t approaches the need to impose this property in 
some arbitrary way was a major source of uncertainty [^,^,^ . 

Another interesting feature of the resummed anomalous dimension is that, for small 
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as, the divergence at uJc is proportional to a1 and not to a^. 
1 dojc 



7eff 



uo — uOc dt 



+ 0(1) 

2^,2 ^, li\ 1/3 



=^,2/3 



UJ — ijJc 



a; — cUf. <^ 1 



(5.2a) 
(5.2b) 



which follows from the linear behavior of the regular solution close to the zero, e.g., in 
the Airy representation of Eq. ( [4.7|) . The singularity ( ^.2|) causes the effective splitting 
function to be power behaved for x — 0, i.e., Pefr(a;,t) ~ x"'^'^'-*''. Note however that, 
since co'c(^) comes from a zero of guj{t), the singularity (|5.2|) does not necessarily transfer 
to guj(t) itself which, according to Eq. ( p.28| ), is expected to have an essential singularity 
a.t UJ = only, even if a complete analysis of possible singularities in the complex cj-plane 
is still needed. 
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Figure 6: lOc and lOg as a function of for the BFKL kernel with nj = 0. 



The values of the exponents uJc and uJs as a function of as (and Q"^), are shown in Fig. ^ 
and compared with the LL and pure NLL results. It is apparent that the improved equa- 
tion provides sensible predictions even for sizeable values of Og. A significant difference 
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between the two resummed exponents ujc and ujg persists even to low values of as, largely 
as a consequence of their differing by a slowly convergent series of non-integer powers of 



cts, as discussed in section 4.2. 



The above difference should not be too confusing. The exponent ujs{t) signals the 
breakdown of the formal small-x expansion of the anomalous dimension of Eq. ( p.31| ), 
due to infinite saddle-point fluctuations, while udt) tells us the position of the singular- 
ity of the resummed anomalous dimension. Their difference arises from their different 
definitions, not from some instability of our approach (cf. Sec. 



What is the relation that such quantities bear to the pomeron singularity tup, the 
leading cj-plane singularity of the gluon Green's function? Though the latter is dependent 
on the strong coupling region, we expect that > Maxci;c(^) for a positive definite as(t), 
due to the very definition of ujc(t) as a zero of the integrated regular solution Qujit), to 
which J-'uj{k) is closely related (Sec. W^. In fact tup is defined as the value of u being itself 



equal to the endpoint of the spectrum: up = /ip(c<jp) (Sec. |2l3| ), and thus corresponds to 
a nodeless jF^(fc), regular for t — > — oo also. Therefore, if the interaction does not change 
sign {as(t) > 0), J-'ujik) can have a node for < cup only, so that ujcit) < cup. 

The above remark implies that the small-x behavior of the gluon Green's function, 
dominated by the singularity ai uo = uo^ in J^(^(fco), is not sensitive to the region uj ~ ujf.{t) 
where J^ui{k) changes sign. This fact is consistent with the positivity constraint on the 
total cross section. 



maximizes 



Furthermore, we can state that the frozen regularization of Sec. ^[3 
the interaction strength in the strong coupling region t compared to various cut-off 
procedures. Therefore we also expect ujf < ojgit), the value quoted in Eq. (|2.18| ) at the 



freezing point. It follows that Maxti;c(t) ^ ^^p ^ Maxco's(t) or, in other words, that the 
two exponents of Fig. ^ provide, in the strong coupling region, lower and upper bounds 
on the pomeron intercept cup. Of course, the precise value of the latter will be dependent 
on the size and shape of the effective coupling in the small-fc region. 
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5.2 Stability 



The original LL + NLL formalism suffered from considerable instabilities under renormal- 
ization group scale and scheme changes. 

An important characteristic of any resummed approach is that it should be relatively 
insensitive to such changes, and generally stable. In the approach advocated here, it has 
already been shown in the previous sections that the formal truncation error is small. It 
still remains to demonstrate its stability in practice. 



Re normalization scale and scheme. Note first that in our approach the renormal- 
ization scale only enters through the RG invariant A parameter (Eqs. ( p.2|) and (|2.28| )). 



It is then easy to see that the physical results are A- independent. A redefinition of A is 
essentially a shift in t, say by an amount At. There is a corresponding modification of 
Xi , X2 5 • • • by the amounts 

X-^Xi+bAtXo, X2 +26Atxr + &'(At)\^, ... (5.3) 

In the off-shell 7-representation ([4.2[), this corresponds to a modification of by an 



amount bfi'yAt. In fact the transformation (|5.3|) changes the coefficient rji only, the re- 
maining ones ri2, rj^, . . . being left invariant. This change exactly cancels the modification 
of t itself: 

exp |7t - i-X^(7,/i)| ^ exp !^j{t + At) - l-{X^{j,fi) + bfi-fAt)^ , (5.4) 

thus implying that the physical results are independent of the A-parameter choice. This 
automatic resummation of the renormalization scale alleviates the need for techniques 



such as BLM resummation [Q, advocated for example in ||3^, which show a strong 
renormalization scheme dependence. 

The issue of renormalization scheme dependence is in fact closely related. Consider a 



scheme S related to the MS scheme by 

af ) = ap) + Ta! , (5.5) 
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with an appropriate modification of Xi- Except for terms of O {aD and higher, this is 
identical to a renormahzation scale change. Indeed if one defines the scheme change by a 
modification of A then renormahzation scheme changes behave exactly as renormahzation 
scale changes, and so have no effect on the answer. Using instead ( |5.5| ) there is some 
residual dependence on the scheme at 0{al), but as one can see in Fig. |^ for the T 
scheme, which has T ~ 1.17 (for = 0), the effect of the change of scheme is small. 
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Figure 7: Renormahzation scheme uncertainty of the two exponents; MS scheme and T 
scheme; is always shown in the MS scheme, and is connected to the T scheme value 
of via (|5.5|) . 




Resummation scheme. In resumming the double transverse logarithms (energy-scale 
terms), there is some freedom in one's choice of how to shift the poles around 7 = and 
7 = 1. In a similar manner to what was done in 0] we consider two choices. The one 
explicitly discussed in this paper (and the one used for all the figures elsewhere in this 
paper) can be summarized as 

(5.6) 
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with an equivalent procedure around 7 = 1. We refer to this as resummation type (a). 
An alternative possibility is 

(5.7) 



7" (7 + ^ujy 

Thus we have 



Xo(7)=Xo(7)---^+ ^ 



7 I-7 7 + ^t^ 1-7 + ^u;' 
Xi (7) = Xi(7) + . \ \o + ^ ^ . + - 



(7 + icj)2 (l-^+icu)2 2V7 + i^ 1-7 + ^^ 



(5.8a) 
(5.8b) 



Xi(7) = Xi(7) + 



Xo(7) 



1 1 

+ 



7' (1-7)' 



1 /I 1 

+ 



2 V7 1-7 



Ai(0) Ai(0)-6 
7' (1-7)2 • 

(5.8c) 



A comparison of these two resummation schemes is given in Fig. ^ and the difference 
between them is again reasonably small. 
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Figure 8: Resummation scheme uncertainty of the two exponents. 



Aside from the explicit renormalization-scale independence, the stability of our ap- 
proach is connected with the resummation of the coUinear poles, for both the double-log, 
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energy-scale dependent terms (the 1/7'^ and 1/(1 — 7)'^ poles at NLO) and for the single- 
log ones of Eqs. ( [4.11| ) and ([4.12| ). Stability has been noted elsewhere, in the study of a 
rapidity veto (initially examined in |^) combined with a resummation of the energy-scale 
terms 



6 Conclusions 

In this paper, we have improved the small-x equation in several ways. Firstly, we have 
taken into account the collinear limits, and their scale dependence. This implies the 
cu-shifts of the 7-singularities in Eq. ( |2.6|) , which yield a double- log resummation of pa- 
rameters like C(j/7 ~ as/7^ or 07/1 — 7, and implies also the effective characteristic function 
in Eq. ( |4.12| ), which yields a single- log resummation in the parameter OsXi/Xo — ots/l or 
tts/l - 7- 

Both kinds of resummation require an infinite number of subleading terms in the 
original BFKL formalism, in which both uo and the running coupling play the role of 
expansion parameters. The RG improved kernel (Eq. ( |2.2| )) is actually an infinite series in 
as{t) of cu-dependent kernels, so that the corresponding Eq. (|2.8|) is no longer an evolution 
equation in logl/x with a simple dependence on the conjugate variable a;, but a much 
more general cj-dependent integral equation. 

The second important improvement concerns the treatment of this generalized equa- 
tion. In the limit in which the Green's function is factorized (Eq. ( |2.7| )) we have singled 
out the solutions of the homogeneous equation J-'uj{k) (J^a;(fco)) which are regular for 
if: — > +00 {t —00), and we have provided a general method for the construction of 
Tu}{k) in Eq. ( |2.19| ). The latter exploits u as expansion parameter (or the eigenvalue 
H ^ oj ii referred to the eigenfunctions ( p.ll| )), and thus we call it cu-expansion. It allows 
the construction, in terms of the improved kernels, of the characteristic function x{l^^) 
of Eqs. ( |4.2| ) and ( |4.4| ), which shows no sign of instability when uj increases (Fig. |^), even 



if the improved kernel is truncated at NL accuracy. 
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The key advantages of the improved equation concern the calculation of the t-depen- 
dence, or of the resummed anomalous dimensions, which can be given in terms of T^jik) 
only. Let us list some of them: 

• The resummation involves not only all powers of ag/uj, but also an infinite number 
of subleading terms and extrapolates quite smoothly the fixed order perturbative 
result (Fig. ^). 

• Although we resum only a fraction of such subleading terms, we have characterized 
the error that we make as a constant scale change At = 0{uj), or Aos/Qs = 0{as uj). 
Therefore, the neglected terms are subleading, order by order, in both as/uj and as 
expansions. 

• Although we are limited in principle to small cu's, we incorporate exact one-loop (and 
partly two-loop) anomalous dimensions in the tu-dependent kernels. In particular, 
we have exact energy- momentum conservation, i.e., the gluon anomalous dimension 
vanishes for = 1. 

We have provided two critical exponents ujs{t) and ujc{t) that signal the breakdown of the 
above resummation. The first one {ujs{t)) is roughly related to the breakdown of the ag/uj 
resummation, or better of the saddle point ( "semiclassical" ) approximation, valid for large 
boot (or ag/uj <^ 1), and was the only one considered in previous L-I-NL estimates. The 
latter exponent {u!c{t)) comes from a zero of the gluon density and signals a singularity of 
the resummed anomalous dimension series. Their difference involves non-integer powers 
of as {as and higher) which are related to a "quantum" wavelength in the t-dependence. 

The estimates of ujs{t) and ujc{t) (Fig. P) in the improved formulation are now quite 
stable (Figs. ^ and P) — despite the large size of NL corrections — and nearly renormaliza- 
tion-scheme independent! The reason for that stems from both the coUinear improvement 
of the kernel, and from the RG invariant formulation of the solution. Both exponents 
are actually useful for a full understanding of the solution J^(^(A;o), carrying the (non- 
perturbative) pomeron singularity cjp. Indeed we have argued that — for reasonable strong 
coupling extrapolations (positive definite as(t)) — the pomeron intercept is bounded 
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between Maxoodt) and Ma.xuJs{t). Present estimates of the latter (Fig. ^ are consistent 
with the small-x exponent ^ 0.2 seen for moderate at HERA [Q. But a detailed 



analysis, including two-scale processes [41,42], is required to obtain a clearcut picture. 



Having no problems with stability, we are now more confident of future progress. We 
have already mentioned the need for evaluating jF(^(fc), the regular solution for t — oo, 
which is much more dependent on the strong coupling region. But also the full Green's 
function Quj{k,ko) for k/ko = 0{1) — i.e., outside the the factorization regime — is 



interesting for the description of two-scale processes (double DIS forward jet ||42|| , 
etc.). We hope to have a better understanding of both quantities from a simple model 
with coUinear resummation 



Of course, a complete understanding involves a variety of other questions, like a real- 



istic evaluation of tup, a full inclusion of quarks (Sec. iA) and impact factors [22,43|, the 
relation to the CCFM equation |]5|,0 and other two-channel formulations 0], and so on. 
But we think that, despite some residual uncertainties, we are on the right track. 
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A /i-Expansion 

A.l Saddle-point method 

In this appendix we want to show the saddle point procedure for deriving the dependence 
of the coefficients r^^ : j = 0, 1, . . . of the small- /i expansion for Xuiil: A*) 

xUl, = Voil) + ^2 (7) + ■ ■ ■ , (A.l) 

in terms of the eigenvalue functions : J = 0, 1, . . . of the coefficient kernels in Eq. ( p.. 41 ). 
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The action of the improved kernel on its eigenfunctions is 



= [/C^ - fi]mk) 



27ri 



,n=0 



(A.2) 



Now we assume the above integral to be dominated by a saddle point at 7 = 7j^(/i, t) 
where h^t = x^(7,yu), Xui being the 7-derivative of (see Eq. (|2.20|) ). By adopting 7 
and /i as independent variables, we replace 



as{t) = - 



Introducing the "mean value" 

/ v4(7)e-^"(^''^) d7 



(A.3) 



-K.(7,/^) 



00 

[7^7]^-^E;;^xL™-H7,.)A^ 



(A.4) 
(A.5) 



m=2 



where xL™'' = dll^Xw and A = 7 — 7, we write Eq. ( |A.2| ) in the form 



X(7,/^) 



(Xn(7)) 



n=0 

having dropped the dependence. We observe that 
= 



(A.6) 



lim — lim ) e ^ 

7— »+ioo 7^— ioo 



d7 d, (e-^) oc {d,V) = -[(x(7,/i)) - bfit] 



Collecting Eqs. ( |A.2| ), ( [A.6| ) and ( |A.7| ) we obtain the basic equation 



(X(7,/^)) = E'"'(^^-(^)) = E 



n=0 



.X(7,/^) 



(Xn(7)) 



(A.7) 



(A.8) 



At lowest order we have 



(^0(7)) = (Xo(7)) • 
One can easily check thatQ (A") = Since 



(A.9) 



(^(7)) = (^(7)) + (^'(7)A) + 
^We denote as [x] the integer part of x. 



= A(7) + 0(/i) (A.IO) 
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it follows that, for all 7, rjoi^j) = Xoil) and hence r/o = Xo- Taking into account Eq. ( |A.9D 
we can simplify Eq. ( [A.8| ) as 

,n-l 



iXnb)). 



The lowest order of this new relation yields 

(.hW> = ^ = ^ + o(,) = ^ + o(,) 

X(7,/^) ^0(7) Xo(7) 
and hence 771 = Xi/Xo- The next order reads 

Wi(7)) + /^W2(7)) = 77— A + ■ 



X(7,/^) [x(7,^)]^' 
By expanding with respect to 7 around 7 yields 



(A.ll) 



(A.12) 



(A. 13) 



^1(7) + ^1(7) (A) + -r/i'(7)(A2) + /i 7/2(7) 



Xi(7)+x'i(7)(A) + -x'i'(7)(A^) 



X(7,/^) 

To the relevant order in /i, we have 
¥ Xo(7) 



+ /i 



X2(7) 
[x(7,/i)]^ 



(A. 14) 



(A) 



2 [x'o(7)]^ 



6/i 



X(7,/^) Xo(7) 



1 ( . Viil) 
Xo(7) 



(A.15) 



and substituting in Eq. ( |A.14| ) we get 



X2 , Vt 

V2- — + — 

Xo Xo 



2 [x',f 2x'o \ Xl 



xl 



xl 



— Til 

Xo 



I.e., 



V2 



1 

Xo 



X2 

Xo 



- (?7i - bd^)r]i 



(A.16) 



Going further requires taking into account higher order terms both in the fluctuations 
(A™) and in the /x-expansion of Eq. ( [A. 111 ). 

The advantage of this method is that it is clearly local in t (because the saddle point 
% is a function of t) and in 7 (because of the finite fluctuations). Therefore, if t is 
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large enough for a stable saddle point to exist, then the procedure and the result are 
independent of the regularization procedure in the strong coupling region t ~ 0. 

The disadvantage, though, is that the order of fluctuations required increases rapidly 
with the yU-exponent. It turns out in fact that, in order to determine rjj : j > 2 — i.e., to 
evaluate Eq. (|A.11|) to order — , the most involved calculation concerns (A) which 



requires the computation of the fluctuations in / Ae ^d7 up to order 6j — 8. 
A. 2 7-Derivative method 

By comparison, the 7-derivative method is formally much simpler. We start rewriting the 
basic equation ( p.24| ) by introducing the notation 

D = {r]^-bd,) (A.17) 
in the form 

^0.(7, y") = [Xo + f^D]-' xi + [Xo + /iI^]~V X2 + ■ ■ • 

= f 1 + ^ ^ V74 + Xo-' (i+^Dy\4+...] . (A.18) 

V X.. / \Xo V X<^ / Xo / 

We then expand in fi both the Xn series and the operator denominators, which depend 
on ?7tj in a non-linear way, and we derive the /i-expansion (|2.25|) . 



For instance, if we want r/g , we can rewrite Eq. ( A.18| ) up to order /x^ in the form 



V = (1 -n'D + Axo'D)') ^ + - (1 - /x(Dxo 1 + Xo '^)) ^- + /i'- 

Xo Xo Xo Xo 

D = Di+ fxr]2 + -- - , Di = rii-bd^, (A.19) 

where the u index has been dropped. 

We then identify the rji coefficients in Eq. ( |A.19D term by term: 



Xi ^ fX2 Xi\ /A on^ 

Vi = V2 = — [ Di— , (A.20) 

Xo Xo VXo Xo/ 

X3 1 /^n ^ , ^ n ^ ^2 , / 1 n V Xi 1 Xi 
V3 = Di \ Di h —Di r]2— , 

Xo Xo V Xo Xo / Xo VXo / Xo Xo Xo 
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Eq. i WM ) proves Eq. {^2^ of the text 



We notice the curious fact that if 



X2 ^ Xl X3 r^2Xl 



^ = D^^, ^ = Df^, (A.21) 
Xo Xo Xo Xo 

both ri2 and r^s vanish identically. This is a particular case of Eq. ( |4.1CI| ), which states that 
rji = Xi/Xo is an exact solution of Eq. (|A. 18|) if 

^ = D"^, n>l. (A.22) 
Xo Xo 

In fact we have the chain of identities 

oo 

Xo ^ ^ ' Xo 



^(Xo + fiD.r^-^'hoif^D^r- = J2^Xo + /i^i)-^"+'^Xn+i , (A.24) 



oo 



n=0 n=0 



which prove Eq. ( [7081) if Eq. ( |0^ ) is satisfied. 



It is straightforward to check that the ansatz (|A.22|) builds up the correct collinear 
singularities to all orders. We start from 

, (7-0,1) (A.25) 

Xo 7 1-7 

and, by applying the Di operator of Eq. ( |A.22| ) in the relevant limits we obtain the result 



Xo V7 /7 vl-7 /1-7 

which checks with Eq. (|2.6| ) . It follows that the leading collinear singularities must cancel 



out in r]j : j > 2, as stated in Sec. ^ 



We should keep in mind that the two methods just illustrated are equivalent when both 
make sense, i.e., for t large enough for the stable saddle point to exist. This assumption is 
implicitly present in the 7-derivative method when we expand the operators in Eq. ([A. 181). 



This means that we stay away from the zero modes of the full operator and we consider 
the D operator as a small perturbation with respect to xo- Expanding in D is analogous 
to the fluctuation expansion. 
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